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The transverse anisotropy constant and the related Doring mass density are key parameters of 
the one-dimensional model to describe the motion of magnetic domain walls. So far, no general 
framework is available to determine these quantities from static characterizations such as magne- 
tometry measurements. Here, we derive a universal analytical expression to calculate the transverse 
anisotropy constant for the important class of perpendicular magnetic multilayers. All the required 
input parameters of the model, such as the number of repeats, the thickness of a single magnetic 
layer, and the layer periodicity, as well as the effective perpendicular anisotropy, the saturation mag¬ 
netization, and the static domain wall width are accessible by static sample characterizations. We 
apply our model to a widely used multilayer system and find that the effective transverse anisotropy 
constant is a factor 7 different from the when using the conventional approximations, showing the 
importance of using our analysis scheme. 

PACS numbers: 75.78.Fg,75.70.-i 


I. INTRODUCTION 


Multilayers with perpendicular magnetic anisotropy 
(PMA) are widely used in research as well as in 
applications. 1 A particular focus of the research today 
is the investigation of domain wall dynamics in such ma¬ 
terials. Theoretically, the dynamics of domain walls is 
often described 2 by the one-dimensional model. 3 In this 
model, the tilt angle ip of the spins in the domain wall 
with respect to the domain wall plane is used as the con¬ 
jugated momentum to the domain wall position q. The 
energy associated with a change of ip can be described by 
an effective uniaxial anisotropy, the so-called transverse 
anisotropy with anisotropy constant K±. In general, this 
anisotropy acts as an energy reservoir that leads to do¬ 
main wall quasiparticle behavior, such as domain wall 
inertia 4 due to an effective domain wall mass. 5 

The transverse anisotropy constant is a key parame¬ 
ter to describe domain wall dynamics: it is proportional 
the to critical excitation strength (field or current) be¬ 
yond which the domain wall becomes non-steady (Walker 
breakdown) 2 and it is related to the inertia of the domain 
wall through the Doring mass density m o through the 
simple relation 3,6 


m D 


1 + a 2 ) 
K± 7 2 Aq ’ 


(1) 


where a is the Gilbert damping, 7 = 1.76 x 10 11 As/kg 
is the gyromagnetic ratio, and Ao the static domain wall 
width. 

The effective transverse anisotropy constant K± is 
not directly accessible from static sample characteriza¬ 
tion. Existing theories predict K± only for homogeneous 
magnetic materials and only in the limits of infinitely 


thick samples (bulk samples), 3 where I\± = 1/2/zoM 2 , 
and for ultra-thin films of thickness T —> 0, 7 where 
K± = ln(2)T/xoM 2 /(27rAo). As shown in this paper, 
these approximations yield significantly inaccurate re¬ 
sults if applied to typical multilayer films, which typically 
are the application-relevant systems. 

In this paper, we present rigorous calculations of the 
magnetostatic energy of a multilayer system with a do¬ 
main wall, from which we obtain analytical expressions 
for the transverse anisotropy constant K± as a function 
of the thickness of a single magnetic layer T, the multi¬ 
layer periodicity V, the number of repeats TV, the static 
domain wall width Ao, and the saturation magnetization 
M s . We find that the saturation magnetization enters the 
equation only as a simple linear scaling factor and that 
Ao, V, and T enter only in the ratios p := V/(2ttAo), 
t := T/(27tA 0 ), and r := V/T. We provide accurate 
equations for K± as the main result of this paper in 
Eq. (40) and we will show that 


K±(t,p,Af) 


PqM 2 


d 2 f(r,Af) + dln(2) 
2d 2 f(r,M) +dg(r, Af) + 1 


( 2 ) 


with 


/(r, AO = (0.9r — 0.76) /J\T 

g(r,N) = r(1.8 - 0.05 ln(7V + 12) + 0.05 ln(r)) (3) 
d = tAf 


approximates the exact result with less than 9 % error for 
10“ 4 <W< 10 5 , 1 < A f < 128, and 1.25 < r < 8 and 
for r = 1 and N = 1, which covers many systems that 
are intensely investigated today. 
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i = 0 

Figure 1. Magnetization profile assumed by the calculations 
in this paper. The profile shows a domain wall at x = 0, sep¬ 
arating a domain where all magnetic moments point upwards 
(black arrows at the left side of the picture) from a domain 
where all magnetic moments point downwards (white arrows 
at the right side of the picture). The angle in the z-y-plane 
is the so-called transverse angle ip. This angle is constant for 
all magnetic moments. 


derivative 1 . We consider a thin film multilayer sys¬ 
tem of length C (cc-direction) and width W (y-direction) 
in the limit £, W —> oo, which is a reasonably accurate 
approximation if C, W Ao- Temporarily, we will treat 
for the calculation the width as a finite value and set it 
to infinity at a later stage. The ground state of such a 
perpendicular magnetic film is characterized by s = 0. 
As in the ID model, we will assume the presence of a 
single domain wall in the y-z-plane at x = 0. However, 
now we have the additional y dependency of the finite 
width and the z-dependency of the multilayer, which we 
describe as follows: 


II. GENERAL EQUATIONS 


The transverse anisotropy constant is a parameters of 
the one dimensional (ID) model. The following calcula¬ 
tions are based on this model, which assumes the pres¬ 
ence of a single, straight domain wall in the y-z-plane 
with the magnetization profile 

M x (x) = M s cosh _1 (a:/A) sin(i/’), (4) 

M y (x) = M s cosh _1 (a:/A) cos(ip), (5) 

M z (x) = M s tanh(x/A), (6) 

as depicted in Fig. 1. Here, M = (M x , M y , M z ) is the 
local magnetization, M s = |M| is the saturation magne¬ 
tization, A = A (ip) is the domain wall width parameter, 
and ip is the transverse (domain wall) angle. The magne¬ 
tization is homogeneous along the wall, i.e., it does not 
depend on y and z. The magnetostatic energy density a d 
of such a domain wall (energy per unit area of the wall) 
is symmetric in ip and therefore can be described as a 
function of s := sin 2 ip. Also, the magnetostatic energy 
can be read as a uniaxial anisotropy with the hard axis 
along the x direction. The general formula for such a 
uniaxial anisotropy reads 

a d {s) =C + I<± I d x{M x /M s ) 2 + 0(s 2 ) (7) 

= C + 2A 0 K±s + 0(s 2 ), (8) 


where C is a constant and K± is the leading order 
anisotropy constant. The first order partial derivative 
of Eq. (8) at s = 0 yields K±: 


K± = 


& 

2A 0 ’ 


( 9 ) 


where the superscript < - 1 - ) denotes the partial derivative 
with respect to s at s = 0. We will hence derive a formula 
to calculate the magnetostatic energy 


M-l 


M x (x, z) = M x {x)w{y) ^2 v ( z ~ 

3 =o 
A/"—1 

My(x, z ) = M y (x)w(y) '^2 v ( z - 

3=0 

AT-l 

M z (x,z) = M z (x)w(y) ^2 v ( z ~ 

3=0 

w(y) = 9(W/2-\y-W/2\), 
v(z) = 0(T/2-\z-T/2\). 


jV), 

(10) 

jV), 

(11) 

jV), 

(12) 


(13) 


(14) 


The magnetostatic energy of a magnetization configura¬ 
tion M(r) per unit domain wall area reads 


aiis) = tew(vT // d 3 rd 3 rV(r) / ,(r') F A_, (15 ) 


where p is the volume density of magnetic charges, W 
is the width of the structure (y-direction) and AfT is 
the total thickness of the magnetic material (z-direction). 
Magnetic charges arise from the divergence of the mag¬ 
netization inside the volume (volume charges) and from 
the components of the magnetization perpendicular to 
the surfaces of the structure (surface charges). However, 
the surface charges do not explicitly depend on s and the 
interactions between volume charges and surface charges 
vanish when averaged over y and z for symmetry reasons 
since all charges are antisymmetric around r = 0. There¬ 
fore, the only relevant charges are the volume charges 


M-i 

P( r ) = ^2 P^x,y,z - jV), 

3=0 


pi(r) = Pl {x)w(y)t(z), 


pi (x) = — div M = sin('0) 


tanh(x/A) 
cosh(a;/A) ’ 


(16) 

(17) 

(18) 


III. MAGNETOSTATIC ENERGY DENSITY 

In the following, we calculate the magnetostatic en¬ 
ergy density a d and, in particular, the first order partial 


To solve the stray field integral Eq. (15), we will first sim¬ 
plify the general form of the integral using the symmetry 
properties of the charge distribution Eq. (18) and use the 
explicit form of the charge distribution at a later stage. 
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The kernel of the double sum in the integral 
Mo 


<?d = 


d 3 rdV 


8ttWA/T . 

N-i 1 

X E Pl{x,y,Z - lP)px(x',,z' - jV)~ - yv 

z —< r — r' 

i,j —0 1 1 


Mo 

" 8ttWA fT 

M-l 

><E 


d 3 rd 3 r'pi (a;, y, z)p ± (x ', y', z') 


1 


( 20 ) 


i,j=0 


r + i?e 2 — r' — j”Pe~ 


depends on i and j only in the form i — j. That is, we 
can reorder the double sum to just one over this difference 
(i — j) = (—A f + 1), Af — 1 (which we will index again 
by j) and a factor Af—\j\ that counts how often the term 
j is represented in the original sum: 


<?d = 


Mo 

8tt>VA/T 

M-l 

x E 

j=-M +1 


JJ d 3 rd 3 r'pi(r)pi (r') 


AA-|i| 


( 21 ) 


r — r' + jVe z 


The charge density p has only a trivial dependence on 
y and z. Hence, we can easily solve the integration with 
respect to y and y', giving 


<Jd = 


Mo 


JJ dxdx J dz J dz p\{x)p\{x) 


8tta rr. 

AT-l 

x Y (A - - lil) /w( X - x', z - z' + jT) 

i=-AT+i 


( 22 ) 


with 


/■w 

r w 

/ d y 

/ d y' 

Jo 

Jo 


\/(y — y ') 2 + £ 2 + ~ 2 


= -2 + 2 ln(2W) - In (a; 2 + z 2 ) + 0(W~ 


(23) 


(24) 


The steps to solve the integral Eq. (23) are discussed 
in detail in appendix A. The terms constant in x vanish 
with the integration over x because the magnetic charges 
are antisymmetric in x. That is, in the limit W — > oo, 


= 


Mo 


8ttA/T 


J\f 1 /-T r'T r r 

(A/ - — | j |) J dz J dz' J J dxdx' 


j=~M+ 1 Jo 170 

x pi (a:)pi ( x')h{x — x',z — z' + jV ) 


with 


h(x, z) = — ln(x 2 + 2 2 ) 


(25) 

(26) 


As outlined in appendix B, the double integral over x and 
x' can be reduced to a single integral in Fourier space. 


The result reads (the hat denotes a Fourier transform 
with respect to the ^-coordinate): 


<“> E>-i wfafa 




(27) 


x / dk\s(k)\ 2 h(k, z - z' + jV), 


with 


1 tanh(x/A 0 ) 
A 0 cosh(a;/Ao) ! 


?(a:) = 


?(&) = ik A 0 W - 


2 cosh(^)’ 


h(k,z) = \PhxJ—e 

|fc| 


(28) 

(29) 

(30) 


The only term that depends on 2 is the function h. The 
integration of this part with respect to 2 and z' is lengthy 
but not very involved, as shown in detail in appendix C. 
The result reads 


( 1 ) Mo-M^Ao ^ M - \j\ I., 

^ = —-— Y \r E - 2 ) 

j=—J\f +1 i=— 1 


St 


f°° e -\jp+it\q + \j p + it \ q - 1 

Jo q gcosh 2 (g/4) 


(31) 


where we have introduced the reduced variables q := 

2ttA 0 fc, p := 2^7’ and * : = 2 tE‘ The _1 in the in " 
tegral has been added to make each integral finite; it 
cancels out in the sum over i. We have furthermore re¬ 
placed the integral over all q by twice the integral over q 
from 0 to 00 . The integral 


/>00 

G(x) := / dq 
Jo 


e xq + xq — 1 
0 ’ q cosh 2 (g/4) 


(32) 


can be solved analytically, as shown in appendix D. The 
final result reads 


(l) _ MoMs 2 Ao 

a d — — 


M-\j\ 


8t E tv" 

j=-N +1 


(33) 


x (G(\jp + t|) + G(\jp - t\) - 2G(\jp\)), 


where the function G(x) is explicitly given in Eq. (DIO). 
This result is exact but very lengthy and difficult to eval¬ 
uate. We will optimize the result for numerical evaluation 
in the next section. 


IV. EVALUATION 

For practical purposes, the evaluation of Eq. (33) is 
often computationally too expensive. Also, we note that 
Eq. (33) contains differences of G, which means that a 
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numerical evaluation of G may yield dramatically wrong 
results just due to the finite precision of computer algebra 
and the error made cannot be estimated easily. We will 
therefore re-write Eq. (33) to eliminate the differences 
to make it more robust for numerical evaluation. In the 
following steps, we use that jp > t for j >0 and that 
G(0) = 0 to re-write the sum over j and to eliminate the 
absolute signs in the arguments of G. Subsequently, we 
expand G(x) in a Taylor series around jp. 


= 2 G{t) 
AC-1 


(G{\jp + 1|) + G(\jp - 1 1) - 2G(\jp\)) 

(34) 


+ 2 


3=1 


Af - j 
Af 


(' G{jp + t) + G(jp -t)- 2 G(jp)) 


= 2G(t) + 4 


A/”—1 . r 

\ ' Af-j 


Af 


E 


t 


2m 


(2m)! 


3 =1 

U ~ l \f n °° t 2m 

= 2GW+4VV Et^TT 

^ Af ^ (2m)! 

j —1 m =1 ' 


G {2m \jp) (35) 


F (2m - 1] {jp), (36) 


where we have introduced the derivative of G: 

1 - e~ xq 


F(x) := f 
Jo 


d q- 


= 8a; 


(37) 

(38) 


‘cosh 2 (g/4) 

{«■ - . 

where H x is the harmonic number (see appendix D for the 
steps to solve the integral). The Taylor series converges 
for all t, p, and M and the derivatives F^ ( x ) are easy to 
calculate and to evaluate. It is useful to note that F^ ( x ) 
are all positive (negative) for odd (even) n , finite for x £ 
[0,oo], and monotonically approaching 0 as x —> oo. We 
estimate the error by the Lagrange error bound for Taylor 
series, which states that the error E n by approximating 
a function / at a position x using a Taylor series of order 
n around Xq has an upper bound of 

max |/ (n+1) (:r)| 

£„ < I6K ( I ) + 1) , -1* - (39) 

Inserting this back, we obtain the formula for approxi¬ 
mating up to order M 

IM)M 2 A 0 


r (l) 


(M) = 


41 


(40) 


G(t) 


with the error 

Zb 


At—1 

2E 

3 =1 


A f - j 

u 


M - 1 . 2 m 

T — 

^ (2m) 
m=l v } 


F (2m-l) Up) 


Aa^(M) = 


PoMjAo 


2 1 


(41) 


At—1 

><E 

3 = 1 


M j t2M+1 I f( 2M ) (jp - t ) I 

A r (2M + 1)! 1 [JP )l 



tAf 

Figure 2. Transverse anisotropy constant K± as a function 
of the reduced single layer thickness t, the reduced layer pe¬ 
riodicity p, and the total number of layers Af (points). The 
solid lines show the approximations according to Eq. (2). The 
dashed line shows the approximation of a very thick (bulk) 
sample and the dotted line indicates the results of the thin 
film approximation. Both approximations are significantly in¬ 
accurate in the technologically relevant region of multilayers 
of intermediate thickness. 
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Figure 3. Maximum error made by the approximation of 
Eq. (2), compared to the accurate result of Eq. (40), as a 
function of tAf for Af in the range of 1 to 128 and r in the 
range of 1 to 8. 


where we have used that F < ' 2m ' > is monotonic to derive 
the maximum that enters Eq. (39). 

From Eqs. (9) and (40) we can now derive the trans¬ 
verse anisotropy constant. For the general case of vari¬ 
able Af and arbitrary t and p, the anisotropy constant K± 
is plotted in Fig. 2. We find that this hyper-dimensional 
dependency can be very well approximated by the simple 
rational form given in Eq. (2). The error made by this 
approximation is plotted in Fig. 3 for log 10 (tAf) in the 
range of —4 to 5 in steps of 0.05, log 2 Af in the range of 0 
to 7 in unit steps and r in the range of 1 to 8 in steps of 
0.25. The error is always smaller than 8%. That is, for 
most multilayer systems investigated today, Eq. (2) pro¬ 
vides a very good description of the transverse anisotropy 
constant. 





















5 


V. APPLICATIONS 


ACKNOWLEDGMENTS 


We first note that, for a single homogeneous magnetic 
film, Eq. (40) simplifies to 


K ± (AT=l) = ^L G (t), (42) 

which converges to K± = /iqM 2 /2 for a bulk material 
(t = oo) and to = ln(2)f/zoM 2 for a very thin film 
(f —> 0), as expected. 

The most important application of our calculations 
is the case of a true multilayer system, i.e., for N > 
1. An example is a Pt(2)/[Co68B 3 2(0.4)/Pt(0.7)]3o/ 
Pt(1.3) (thickness in nm) magnetic multilayer with M s = 
1.19(3) x 10 6 A/m and Ao = 11(2) nm, which has been 
used in the investigation of domain dynamics. 1,8 Enter¬ 
ing this into our calculation, we find that the trans¬ 
verse anisotropy constant is given by 0.07(2)/uo M 2 . This 
value is a factor seven smaller than predicted by the sim¬ 
ple and widely used formula K± = hqM 2 /2 and a fac¬ 
tor 1.7 smaller than predicted by the thin film formula 
K± = ln(2 )tnoMg. This means that quantitative predic¬ 
tions of the domain wall dynamics, such as the Walker 
breakdown, 2 are inaccurate when relying on existing for¬ 
mulas, showing that one needs to take into account our 
calculations for a realistic calculation of the resulting do¬ 
main wall dynamics. 
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Appendix A: ^-integration 


Here, we will calculate the integral 

W pW i 

dv l iy 'v(y-yr <A1) 


1 

Wio 


' — y ') 2 + a 

for large W. The integration with respect to y' yields 
1 


r W 


d y'- 


Jo \/(y-y') 2 + a 
= - In (V (y - y 1 ) 2 + a + y-y'^j 
= -f(y- W) + f(y), 
f{y) = In (vV +a + y^j . 


w 

y'= o 


(A2) 

(A3) 

(A4) 

(A5) 


In conclusion, we have derived an accurate analyti¬ 
cal expression for the transverse anisotropy constant in 
perpendicular magnetic multilayer films as used in a ID 
model description. The transverse anisotropy constant 
can be reliably computed using magnetic properties that 
can be ascertained easily for multilayer films from static 
measurements. We have shown that the results can sig¬ 
nificantly deviate from existing and commonly used the¬ 
ories that oversimplify the actual sample configuration 
and are only valid in the very thin or very thick limit but 
not in the intermediate regime, which is however most 
widely used for devices. We have provided a more accu¬ 
rate simplification yielding sufficiently precise results for 
most systems investigated today that are easily acces¬ 
sible without very involved computation. In particular 
we show that for a commonly used multilayer stack our 
result describes the anisotropy well compared to the con¬ 
ventionally used approximations, which are off by up to a 
factor 7. Our results enable precise modeling of domain 
wall dynamics in systems with perpendicular magnetic 
anisotropy using the ID model, which has previously 
failed due to inaccurate assumptions of the transverse 
anisotropy. 


The antiderivative of / reads 

F(y) = yin ^a + y 2 + y) - ^a + y 2 (A6) 

and the total integral is therefore given by 

1 pW pW 2 

wl dy J 0 dy y/iy-y *) 2 + # (A7) 

= (F(W) +F(-W) - 2F(0))/W 

= In (vVW 2 + 1 + l) - In (vVW 2 + 1 - l) ^ 

— 2i/a/yy 2 + 1 + 2^/a/W 2 
= — 2 + 2 ln(2W) — ln(a) + 0(W~ 1 ). (A9) 


Appendix B: Fourier space identities 

Here, we show how to reduce a typical double inte¬ 
gral over a pair interaction kernel to a single integral in 
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Fourier space (assuming real-valued functions): 


JJ dxdx f(x)h{x')g{x — x') 

(Bl) 

= J d xf(x) j d x'g(x — x')h(x') 


= j dxf(x)(h* g)(x) 

(B2) 

= J d xf(x) J d kh(k)g(k)e~ lkx 

(B3) 

= V2tt J d kh(k)g(k) —y= J dxf{x)e~ lkx 

(B4) 

= \p2jt J d kh(k)g(k) ^/(fc)) 

(B5) 


f d z f dz' h(k, z — z' + jV) (C3) 

Jo Jo 

= ^TT r dz [ T dz'e~\ z - z ' +jvm (C4) 

Fl Jo Jo 

= ^ f T dz f T ° V dz'e-\ z '~ zm (C5) 

\k\ Jo J-jv 

= E ( 3 I*I - 2) (e-w™™ + | jV + iT\\k\) . 

1 1 i=— 1 

(C6) 


Appendix D: q-integration 

In this section we provide the exact solution of the 
Fourier space integral Eq. (32): 


Note that we are using the following definition of the 
Fourier transform: 


/>) 

f(x) 


wJ AxI{I)e " 
wJ Aki(k)e ~ ikI 


(B 6 ) 

(B7) 


Appendix C: ^-integration 


In the following we show how to perform the integra¬ 
tion of Eq. (27). We first note that 



(Cl) 


f°° dq e Xq + xq — 1 =: G(x). 
Jo qcosli 2 (q/4) 

To solve this, we note that 

G{x) = [ F(x')dx' 

Jo 

where the integral 

r 00 1 _ p-xq 

F{x) = I dq —^- 


/o cosh (q/4) 


= 4 — 4 


dq e 


— xq 


p «/2 


/(«) 


(1 + e q / 2 ) 2 


nO O 

= 4 / d q(-x)e~ xq 
Jo 

/»oo 

= 8x d qe~ xq 


a'(q) 

-2 


1 + e «/ 2 
e «/ 2 - 1 

( e ?/ 2 - i )( e < z / 2 + 1 ) 


= 8x dq 


oo e —(a=—1/2 )q _ e -xq 


e q — 1 


= 8 x dy 


yX-l/2 _ yX 

i -y 


= 8x (h x - H x _ i) 


(Dl) 

(D2) 

(D3) 

(D4) 

(D5) 

(D 6 ) 

(D7) 

(D 8 ) 

(D9) 


is a continuous second antiderivative of exp(— k\z\). Fur¬ 
thermore, remember that 



= G(c — b) + G(d — 


z) 

a) — G(c 


(C2) 


a) — G(c — b) 


can be solved through integration by parts. Here, H is 
the harmonic number. We therefore obtain 

G(x) = —8 (^j)(- 2 \x + 1) — ^~ 2) ^ 

— s ln(r(x + 1)) + a;In ^ ^ (DIO) 


if 


d 2 G(z) 
dz 2 


f(z). Using these identities, we see that 


where 7 is gamma function and J) is the digamma func¬ 
tion. 
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